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We consider a zero-range process with two species of interacting particles. The steady state phase 
diagram of this model shows a variety of condensate phases in which a single site contains a finite 
fraction of all the particles in the system. Starting from a homogeneous initial distribution, we study 
^jfy the coarsening dynamics in each of these condensate phases, which is expected to follow a scaling 

law. Random walk arguments are used to predict the coarsening exponents in each condensate 
phase. They are shown to depend on the form of the hop rates and on the symmetry of the hopping 
dynamics. The analytic predictions are found to be in good agreement with the results of Monte 
Q\ ' Carlo simulations. 
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I. INTRODUCTION 



Since the first observation of a condensation transition in the homogeneous zero-range process (ZRP) JJ there has 
been a lot of activity to further study this phenomenon on the level of the steady state 0, and on the level of the 
relaxation dynamics 0, Q . When the density of particles exceeds a critical value, the system has been shown to phase 
separate into a homogeneous background and a condensate which contains a finite fraction of all the particles in the 
system. In the steady state the condensate occupies only a single lattice site, and starting with homogeneous initial 
, conditions, the relaxation dynamics exhibit an interesting coarsening phenomenon. 

Besides being of interest in its own right as an example of a condensation transition in an exactly solvable model, 
I | the phenomenon is relevant in a more general context, providing a criterion for phase separation in driven diffusive 
systems 0, || . The basic condensation mechanism is by now well understood on a static and dynamic level, but 
generalizations continue to be a topic of current interest, such as coarsening behaviour on scale- free networks [||, 
processes with defect sites Q or applications to bipartite graphs Of particular interest are generalisations to 
i__T two-species zero-range processes with two conservation laws, which also exhibit condensation and have a much richer 
stationary phase diagram than that of the single species system fil Il0|. Indeed, while the stationary and dynamical 
properties of one-dimensional driven diffusive systems with one species of particles are relatively well understood, 
much less is known about the properties, and in particular the dynamical properties, of driven systems with two or 
■ more species of conserved particles (see [llj for a recent review), 
ly-j ' This paper provides a first analysis of the coarsening dynamics of a two-species zero-range process. We generalize 
the arguments in [2j for a single species system, which turns out to be far from straightforward since several new 
t— I ' effects have to be taken into account, effects due to the coupled dynamics of the two particle species. The model is 
chosen such that all the expected new features can be observed while the steady state is exactly solvable. In addition 
to this theoretical interest, the results are relevant for physical realisations of two species zero-range processes, which 
can be found for example in shaken bidisperse granular systems |12| and models of directed networks |13|. 

In Section II, we define the model which is a generalisation of the model considered in |icj . recap some known 
results for the steady state and give the phase diagram. In Section III we state the expected scaling behaviour for the 
coarsening regime and explain the random walk arguments for its analysis. The main results of the paper are derived 
in Section IV: scaling laws for the time evolution of the mean condensate size for all regions of the phase diagram, 
generalizing the derivation in Q. The predictions are compared to Monte Carlo simulation data and we find good 
O ■ agreement. We conclude in Section V and include a discussion of finite size effects in an appendix. 



II. MODEL 
A. Definition and steady state 



We define the two-species zero-range process on a one-dimensional lattice containing L sites with periodic boundary 
conditions. On this lattice, there are Ni particles of species 1 and N2 particles of species 2. A site with occupation 
numbers k\ and A: 2 for species 1 and 2 respectively, loses a particle of species 1 with rate 51(^1,^2) and of species 2 
with rate g2{k\, k?)- For simplicity we assume that particles hop to their nearest neighbour site to the right, although 
our results also apply for more general hopping of finite range. 

The steady states for this model with general <7i(/ci,/c2) and 32(^1,^2) have been characterised in 0,^3 and we 
now summarise the main points. We denote a particle configuration by k = fci,i, ^2,1! . . . ; Ai^l, hi,L- The steady state 
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probabilities assume a factorised form 

L 

^(k) = n^(k ltX ,k 2 , x ) , (i) 



x=l 

provided the hop rates satisfy the constraint 

5i(fci,fc 2 ) 92(ki,k 2 



9i(ki,k 2 -l) 52(^1-1,^2 
for all fci, fc 2 > 1- The single-site distribution has the form 



(2) 



f.(fci,*a) = ^y/(fci^2)^ fel 4 2 . (3) 
where /(fci, fc 2 ) is a stationary weight which can be written 

/(fci, fc) = fl T^ - ft — 7T — T ■ (4) 
f = i 5i(*,0) ai g 2 (k 1 ,j) 

Here z = (zi, z 2 ), -z, > play the role of fugacities for each species, in that they are chosen to fix the particle densities 
Pi = {Ni) v /L for species i — 1,2, i.e. the expected number of particles per site in the steady state. Thus we are 
working in a grand canonical ensemble, which is normalised by the single site partition function Z(z). This steady 
state can be directly obtained by substitution into the balance condition for the steady state probability that the 
system is in a configuration k. 

We remark that one gets the same steady state if the hopping dynamics are symmetric, rather than asymmetric as 
defined above. A useful property of the steady state is that the expectation value of the hop rate of species i, 

denoted by (gi) u , is equal to z^. Thus {gi) v is a translation invariant quantity; this is obvious in the case of asymmetric 
dynamics where {gi) v is the current, but less obvious in the case of symmetric dynamics with vanishing current, where 

{gi)v > 0. 

We are interested in the coarsening dynamics of the model in various phases that arise for a particular choice of 
rates, namely 

•l + b/(fci + l )T ' 
1 + b/kj 

g 2 (h,k 2 ) = i + b/(h + i)\ (5) 

where gi(0, fc 2 ) = <? 2 (fci, 0) = and 6, c, 7 > 0. It is easy to check by substitution that these rates satisfy the constraint 
(0). The stationary weights, obtained from @, are given by 

fci ' / b \ 

f(h,k 2 ) = — L^(i+ — ) , (6) 



(1+C) fel \ (fcl + l)^ 

where (a)& = Ili=o 1 ( a + *) ^ s the Pochhammer symbol. The single-site partition function is given by 

We make the choice in order to study the behaviour when the dynamics of one of the particle species, here 
species 2, depends only on the number of particles of the other species at the departure site. So condensation of 
species 2, when it occurs, is induced by the presence of species 1 particles, which can be interpreted as an evolving 
disordered background as discussed for a specific case in Section IV. C. The fc 2 dependence in gi is then determined 
by the constraint @. The second factor (1 + c/fci) in gi could be replaced by any function of fci and the steady 
state will still factorise. The form we have chosen is the simplest form of the hop rate for which the single-species 
zero-range process exhibits a condensation transition for c > 2 at a finite critical density of particles Thus the 
parameter c can be tuned to allow also condensation of species 1 particles, which influences the phase diagram of the 
process as discussed in Section II. B. We remark that the existence of condensation transitions and any subsequent 
coarsening behaviour depend only on the asymptotic forms of the rates. Other choices of this second factor, with 
different asymptotic properties, lead either to no transition of species 1 particles if it is non-decreasing or tends to 
a constant faster than 2/fci as fci — > 00, or to condensation at any density (where the fraction of particles in the 
condensate is equal to one) if it tends to zero as fci — > 00. Thus JSJ are basic rates which capture two different 
mechanisms of condensation transition (i.e. induced and autonomous) and the most interesting coarsening behaviour 
that we expect to observe while the steady state remains exactly solvable. 
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FIG. 1: Phase diagram for the choice of rates with b = 1, c = 4 and several values of 7. For 7 = 0.5 and 7 = 1 it is 
c > max(2 + 7, 1 + 27) and region C exists. Whereas for 7 = 2 we have c < max(2 + 7, 1 + 27) and region C does not exist. 
See text for details. 



B. Stationary phase diagram 

The range of possible fugacities is given by the domain of convergence of the partition function Z(z) given in Q. In 
the present case the maximal fugacities are z\ = 1 and Z2 = 1 and when one or both of the fugacities are maximal we 
use the notation z = z c . The phase diagram in terms of the particle densities pi and pi has been derived in |l0lll5| . 
For the grand canonical ensemble the densities are given by 

dlnZ(z) 

Pi = Zi W , i = l,2, (8) 

oz% 

and thus the convergence properties of the partition function at the maximal fugacities determine whether or not the 
critical densities p;. c := Pi\ z ._ l are finite or infinite. In general, pi. c can depend on pi i.e. pi iC = p\. c {pi) and vice 
versa. If pi < pi. c for i = 1, 2, both species are in a fluid phase corresponding to a factorised steady state f z as given 
in J3J. In the phase diagram shown in Figure Q this region is denoted by D. If the particle density pi of either species 
i = l,2 exceeds its critical value p.; jC , species i condenses: the system phase separates into a homogeneous background 
fluid phase with distribution v Za , and a condensate which contains the (p, — Pi, c )L 'excess' particles of species i. In a 
typical stationary configuration this condensate occupies a single, randomly located site. 

Depending on the values of c and 7 in JSJ the following phases appear in the phase diagram in addition to the fluid 
phase D: 

• In region A, the fugacities are z c = (1,2:2) with z^ < 1. The species 1 particles condense and the species 2 
particles form a fluid. The particle densities in the background phase are (pi. c , P"i)- 

• In region B, the fugacities are z c = with Z\ < 1, species 2 condenses and the background particle 
densities are (pi,p2, c )- As an additional point, the site containing the condensate of species 2 particles also 
contains 0(L 1 ^ 1+jS> ) species 1 particles @. 

• In region C, z c = (1, 1). A single site contains condensates of both species and the background densities are 

(Pl,c,P2,c)- 

The phase diagram shown in Figure ^ is richest when c > max(2 + 7, 1 + 27), where all three regions are found. 
For 2 < c < max(2 + 7,1 + 27) the phase diagram contains only the phases A, B and D, and for c < 2 only phases B 
and D remain. 

So far we have discussed straightforward generalisations of previously known results. We now turn to the main aim 
of this work, which is to study the coarsening dynamics of the two-species zero-range process leading to each of the 
condensate phases A, B and C. 



III. COARSENING 



In the following we use the symbol w to denote asymptotic expansions in the thermodynamic limit L — > 00 with 
fixed particle densities, i.e. N± = [piL] and N2 — [P2L]. If the terms in the expansion are only given up to a constant 
factor we use the symbol ~ instead. 
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A. Relaxation dynamics 



In this section we outline the arguments used to describe the coarsening dynamics in the condensate phases. Starting 
from an initially uniform distribution of particles, the dynamics of the condensation can be divided into three regimes: 

(i) nucleation, during which excess particles of either species accumulate at several randomly located sites, which 
we call cluster sites. Each contains O(L) particles, so there are 0(1) cluster sites, separated by a typical 
distance of order L. At the remaining sites, which we call bulk sites, the system relaxes to its steady state 
distribution z/ Zc . 

(ii) coarsening, during which the cluster sites exchange particles through the bulk. This leads to the growth of large 
condensates at the expense of smaller ones and a decrease in the number of cluster sites. 

(iii) saturation, where eventually only two cluster sites remain due to the finite size of the system. In this regime 
the dynamics, under which the system reaches a typical steady state configuration with a single cluster site, is 
different from the coarsening dynamics (cf. 0). 

Physically, the most interesting is the coarsening regime. Here, large condensates gain particles at the expense of 
smaller ones which causes some condensates to disappear. This in turn leads to a decrease in the number of cluster 
sites and hence an increase of the mean condensate size m,i(t), defined as the number of particles of species i = 1, 2 
at cluster sites divided by the number of cluster sites at time t. In the limit t — > oo, rrii(t) converges to its steady 
state value (pi — pi jC )L, the number of excess particles of species i in the system. Within the coarsening regime the 
increase of the mean condensate size is expected to follow a scaling law, (rrii(t)} L ~ t"*. Moreover, on a certain time 
scale r = t(L), the growth of the normalised mean condensate size is expected to be independent of the system size 
L, 

(m l (t)) L ft 



{pi - Pi,c)L 



(t/rr ■ (9) 



Therefore the scale L of the mean condensate size and the time scale r are connected via L ~ t® 1 , and thus r ~ L 1 /^ i . 
The angled brackets denote an ensemble average in a finite system of size L, starting with a homogeneous 

distribution of Nj = [pjL] particles for both species j — 1,2. This is in contrast to the steady state expectation 
denoted by {..) v . The scaling law © defines the exponent which may depend on the particle species i. In general 
one could choose different observables to monitor the coarsening process, such as the square sum of occupation 
numbers. But in our case the mean condensate size is a natural choice, since it is directly accessible by our arguments 
given below. 

We remark that the scaling law is of the same form as that which describes the growth of characteristic length 
scales in phase ordering dynamics [lq . More precisely one can define a scaling function 

(mAt'T)) T 

hi(t) := lmi } " L , (10) 

L^oo (pi — Pi^ c )L 

for all if > 0. With the appropriate time scale t, which will be derived in the next section for the various phases, hi 
is expected to be a non-degenerate, smoothly increasing function with the asymptotic properties 

hi(t') = 0(t' ft ) for t' -> and lim h,-(i') = 1 . (11) 

t' — >oo 

So for small if the coarsening regime is described by a power law J3J which we study in the following. For if — > oo the 
system saturates and hi converges to its maximal value 1. We do not further discuss the behaviour in this regime, 
this has been done for a single species system in 0. 



B. Random walk arguments 

In the following, our aim is to estimate the exponents (3i in each of the condensate phases A, B and C of the model 
defined by the rates This is achieved by adapting the random walk arguments given for the coarsening dynamics 
of the one-species model ;2j. They are based on two major assumptions, which are self-consistent and confirmed by 
simulation data. 
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(Al) Separation of time scales: 

The nucleation process is very fast so that during the coarsening regime the bulk sites have already relaxed to 
the steady state distribution u Zc , 

Within the coarsening regime the system can therefore be separated into a stationary bulk and a finite number of 
isolated cluster sites. On top of stationary hop rates (gi) v = £j )C) cluster sites of species i exchange particles through 
the bulk on a slower time scale, given below. The bulk can be seen as a homogeneous medium through which these 
excess particles perform a biased random walk, and the cluster sites as boundaries where they enter and exit. 

(A2) Independence of excess particles in the bulk: 

The excess particles exchanged by cluster sites perform independent (biased) random walks through the bulk 
on their way to the next cluster site and do not effect the bulk distribution v Zc . 

This is justified below by noting that the average density of excess particles in the bulk vanishes for L — > oo. 

The random walk argument then proceeds as follows. We consider the case where one species i = 1 or 2 condenses. 
The rates we consider decay as <?, — 1 ~ k^ a with < a < I (see Section IV) and the average hop rate in the bulk is 
z i c = 1. So the effective rate at which cluster sites with k t ~ L lose particles to the bulk is gi — z i c ~ L~ a . These 
excess particles perform a biased random walk through the bulk with drift 

(gi i h > o), - ( 9i ) v = - — (g ; } ; - ( 9l )„ = - — l n — -r i . (12) 

Since v Zc (ki — 0) > this is positive and independent of L. Thus the time it takes an excess particle to reach a 
neighbouring cluster site scales as the typical distance between cluster sites which is O(L). So n independent excess 
particles exit the bulk with rate 0(n/L) which has to balance the entry rate of order L~ a . Hence, the number of 
excess particles in the bulk scales as 0(L 1_Q ) which grows only sublinearly with L for < a < 1, consistent with 
(A2). 

In this balanced situation the time scale on which cluster sites exchange single particles through the bulk is 0(L a ). 
The time scale on which cluster sites exchange a finite fraction ~ L of their particles is thus 0(L 1+a ). Since by 
definition the number of cluster sites during coarsening is of order 1, this sets the coarsening time scale r and the 
coarsening exponent /?, in Q to be 

r~L 1+Q , & = — ^- (13) 
1 + a 

With the above considerations we can give an additional motivation for the scaling behaviour We have seen that 
the rate at which cluster sites exchange particles through the bulk depends on their size as k~ a . Thus, in a very rough 
approximation, the time derivative of the average condensate size (TO,(t)), should be proportional to the average 
exchange rate of excess particles, 

K dt ]L ~ (mi(t)) L a ■ (14) 

As the solution we recover the scaling law @ with exponent Pi as given above. This is of course not a strict argument 
and should not be understood as a derivation of the scaling law. 

Compared to the bulk dynamics the coarsening is a very slow process and typical configurations with cluster sites on 
top of a stationary background are quasi-stationary, i.e. within times of order 1 the configurations at cluster sites do 
not change on average. To leading order in L the dynamics on cluster sites have to be compatible with the stationary 
bulk dynamics. By compatibility we mean that the translation invariance of (gi) v implies that it must be the same 
for all sites (both cluster sites and bulk sites) in the system. For two-component systems, this induces consistency 
relations between the occupation numbers k\ and &2 on cluster sites, a fact that will often be used below. If these 
relations are not fulfilled, the configuration is not quasi-stationary in the above sense and changes on time scales of 
order 1 through interaction with the bulk. 

IV. COARSENING SCALING LAWS 
A. Theoretical predictions 



We now apply the arguments described above to the model at hand and explain additional effects due to the 
presence of two species of particles. 
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1. Phase A 

In phase A, pi > p 1)C and only the first species condenses. There are (pi — pi, c )L excess particles of species 1 in the 
system and at the cluster sites k\ = 0((p\ — pi^ c )L), while k 2 remains finite in the limit L — > oo, which is justified 
below by compatibility with the bulk. Hence, at the cluster sites the rates (JjjJ, up to first order in k\, are given by 

9x{kx,k 2 ) m 1 + c/fci , g 2 (h,k 2 ) « 1 + b/kJ . (15) 

Thus the coarsening of the species 1 particles is independent of the second species: the net rate at which particles 
leave a cluster site is g\ — 1 = cjk\. Following the arguments leading to l|13fl in Section III.B the coarsening time 
scale is thus 

r A ~ [(pi - Pi,c)L} 2 , (16) 
and we expect that the normalised mean condensate size grows like 

( {mi{t) \ ~{t/T A )V> i.e. 01 = 1/2. (17) 
(pi - pi,c)L 

This recovers the known coarsening of the condensate in the one-species ZRP where particles hop with rate 1 + elk 

013 

Further, because the jump rates of both species are coupled, the presence of a species 1 condensate influences the 
distribution P of the species 2 particles on the cluster site: Since cluster sites and bulk have to be compatible, g 2 on 
the cluster site has to be equal to the bulk steady state current (g 2 ) v = z 2 < 1, and using l|15|) we have 

(g 2 )„ « (1 + b/kJ)P(k 2 > 0) » P(k 2 > 0) . (18) 

Therefore P(k 2 = 0) » 1 — {g 2 ) v . This is non-zero, but smaller than the expected bulk value, which contains an extra 
positive contribution due to b/kj = 0(1). 

2. Phase B 

In phase B, p 2 > p 2 , c and the second species condenses. The number of particles at a cluster site is k 2 = 0((p 2 — 
p 2 c )Lj. Now, using {J5J|, the hop rate g\ of the first species at a cluster site vanishes in the limit L — > oo if k\ = 0(1). 
But since in the bulk the mean hop rate of the first species is given by its steady state value {g\) v — Z\ € (0, 1) for 
pi > 0, k\ has to be large at cluster sites. Thus, considering k\ large in (JSJ, the hop rate of species 1 particles at 
cluster sites becomes 

gx(kxM) «exp(-6 7 fc 2 /fc 1 1+7 )(l + c/A :i ) , (19) 
which should be consistent with the expected bulk value (gi} u — This compatibility requirement leads to 

k\ +1 w -trfk 2 /lnz! , (20) 
and this relation between k\ and k 2 is dynamically stable since on cluster sites 

dMh, k 2 ) - zi)\ kl+ ^_ b7k2/lnzi « ' (1 + 7 ^ Z1 lQgZl > . (21) 

So cluster sites where ki is too small gain species 1 particles from the bulk (or if k\ is too high they are lost to the 
bulk), and thus any perturbation of the relationship H2()|) is driven towards this stable form on intermediate time 
scales. Hence cluster sites at which (I2U|I is satisfied dominate the coarsening and the hop rate of the second species 
can be written 

g 2 (fci,fc 2 ) = l + 6/(fci + l) 7 ~l+ f 7fc2 - j ■ (22) 

Therefore particles of species 2 escape from a cluster site at a net rate proportional to XjklJ^ 1 ^ and we can repeat 
the arguments given in Section III.B to deduce the coarsening time scale 

1+27 

t b - \(P2 -p 2 .c)L]— . (23) 
The normalised mean condensate size grows like 



(m 2 (t)) L _ui- n(i+ 7 )/(i+2 7 ) . „ a _ 1+7 



(t/T B )WW + *» , i.e. /3 2 = __L. (24) 



(p2~P2,c)L " °' l + 2 7 



7 



3. Phase C 



In phase C, p\ > p\ iC and p 2 > p 2 , c and both species condense. While in phases A and B the relationship between 
the occupation numbers fci and k 2 was fixed by compatibility with the bulk dynamics, in phase C this relationship is 
not uniquely determined. Using the expansion 



we see that for k\ = O(L) any value of fc 2 in the range 0(1) < fc 2 < O(L), and for fc 2 = O(L) any value of k\ in the 
range 0(L 1 ^ 1+1 " > ) < k% < O(L), lead to g\ « g 2 ~ 1 and are compatible with the bulk dynamics. All compatible 
relations between k\ and fc 2 may be observed during the coarsening regime, but the sites with the longest lived relation 
will determine the coarsening timescale. 

Since the leading order of k\ in the hop rate g\ given in l|25() depends on 7, we have to distinguish three cases: 

7<1: In this case the longest lived relation is given by double cluster sites, i.e. sites with k\ ~ k 2 ~ L. Since 
gx — 1 = —b r yk 2 /kj +1 < such cluster sites gain excess species 1 particles from the bulk rather than losing them. 

Double cluster sites are stable compared to cluster sites with other relationships between k\ and fc 2 , in the 
sense that such sites are driven towards ki ~ fc 2 ~ L: For k 2 = 0{L) and < k x < <D(L) one has 

<7i(fci,L) — gi(L,L) ~ ~bjL/kl + ^. Therefore the smaller the value of fci at cluster sites, the greater the rate 
at which species 1 particles are gained from the bulk. Thus ki is driven towards a value O(L). On the other 
hand, for k\ = 0(L) and k 2 < O(L) one has gi(L,k 2 ) — gi(L.L) ~ &7/L 7 > 0, so cluster sites at which only 
ki = O(L) lose species 1 particles to double cluster sites. 

Since on double cluster sites both species exchange particles with the bulk at an effective rate proportional to 
l/fc/, i = 1, 2, the coarsening timescale is given by t(ki) ~ kik] . Thus tq ~ L 1+7 and both species coarsen with 
the same exponent (3i = 1/(1 + 7), i- e - 



7=1: The longest lived sites are again double cluster sites, at which k\ ~ k 2 ~ L. Now the sign of 51 — 1 « —bj/L + c/L 
depends on 6 and c, but all the arguments for 7 < 1 apply in this case also, so we expect that the scaling law 
(EU still holds for 7 = 1. 

7>1: Now the leading order for cluster sites of the first species changes to g\ ~ 1 + cjk\ independent of 7 and fc 2 . 
Thus species 1 coarsens independently of species 2 with the dynamics determined in the same way as phase A, 
therefore j3\ = 1/2. However, the relation between k\ and fc 2 is not stable on cluster sites at which k\ — 0(L) 
for any value of fc 2 , since 52 — 1 ~ b/kj < gi — 1 ~ c/k±. But when k\ is large, fc 2 is driven towards large values 
(since g 2 — 1 is small), thus excess species 2 particles accumulate at sites where fci = O(L). 

On species 2 cluster sites, i.e. sites where fc 2 ~ L, the slowest timescale in the dynamics of the species 2 particles 
is set when the cluster site contains ki = O(L) species 1 particles. However, since the effective exit rates, 
gi — 1 and g% — 1, differ for each species, the coarsening mechanism is more complicated than in previous cases. 
This can be seen as follows. When k\ = O(L), species 2 particles are lost to the bulk with an effective rate 
proportional to C(L~ 7 ). Now, the time it would take for O(L) species 2 particles to escape to the bulk scales 
as 0(L 1+1 ) which is large (since 7 > 1) compared to the timescale 0(L 2 ) over which the species 1 particles 
coarsen. Therefore after a time of order C(L 2 ), the number of species 2 particles at a cluster site is still 0(L) 
but the number of species 1 particles has decreased to its minimum value allowed by continuity, C(L 1 / < - 1+7 '). 
Now the species 2 particles are lost to the bulk in a time of order 0(_L 1 + 7 /( 1 + 7 )) which is fast relative to the 
time of order 0(L 2 ) we have already waited for the species 1 particles to coarsen. Thus the species 2 cluster 
dismantles immediately following the dissolution of the species 1 cluster. Hence both species coarsen on a time 
scale tc ~ L 2 and we expect 



The coarsening of species 2 almost exclusively takes place on vanishing species 1 cluster sites. In this sense the 
coarsening of the species 2 particles is effectively a slave to that of the species 1 particles. Indeed in simulations 
this picture is confirmed, and both species coarsen on the same time scale, but the species 1 particles coarsen 
first (see Figure 0] in the next section). 

The results for each phase are summarized in Table [I] These theoretical predictions are compared with numerical 
results, which are presented in the next subsection, in Figure [21 



gi (fci, - bjk 2 /kl + ~ l + c/ki , # 2 (fci, k 2 ) «ll b/k 



.7 

1 ! 



(25) 




(26) 




(27) 
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Phase 


Coarsening exponents 


A 
B 

C, 7< 1 
C, 7 >1 


01 = 1/2 
01=^2 = 1/2 



TABLE I: Coarsening exponents for asymmetric hopping 
1 

0.9 
0.8 
yS 0.7 
0.6 
0.5 
0.4 

0.5 1 1.5 2 2.5 3 

r 

FIG. 2: Comparison of the theoretical predictions, indicated by lines, with the numerical estimates of the exponents obtained 
in phase A (A), B (o) and C (□). Errors are of the size of the symbols. 



B. Comparison to simulation data 

The theoretical predictions of the previous subsection are compared to Monte Carlo simulations in Figures [21 and 01 
N% = [piL] resp. N2 — [P2L] particles of species 1 resp. 2 are initially distributed on a lattice of size L with uniform 
probability. Cluster sites of species i are defined by the threshold (pi — p^ c )i/40. The proportionality factor has to 
be chosen such that bulk fluctuations are well separated from cluster sites. Results for exponents are not sensitive to 
this choice for the system sizes considered, ranging from L = 256 to 4096, since the fluctuations grow only sublinearly 
in L. With this threshold we measure the mean condensate size rrii(t) and other observables, such as the bulk density 
Pi,buik(i)j of species i as a function of time, scaled with the expected coarsening time scale r. The ensemble average 
is approximated by averaging over 400 sample runs for each system size. 

In the following we discuss the expected behaviour of the observables which is consistent with simulation results, up 
to finite size effects, discussed in the appendix. Plots of the normalized mean condensate size {nii(t)) 1 / (pi — Pi, c )L for 
different system sizes L against the rescaled time t/r, where r is the predicted coarsening time scale, are expected to 
collapse onto a single curve. Within the coarsening regime this curve should be described by the scaling laws derived 
in Section IV. 

During nucleation, more and more particles become trapped in cluster sites, therefore the bulk density pi.bu\k{t) is a 
decreasing function of time, approaching the critical density pi^ c . This is used as a criterion to identify the beginning 
of the coarsening regime. The end of the coarsening regime is reached approximately when (rrii(t)) l = 0.4 (pi — pi tC )L, 
corresponding to an average of 2.5 remaining cluster sites. For later times the data significantly deviate from the 
scaling law and the system saturates, as already explained in (jTU|l . Within the coarsening time regime defined above 
we make a linear fit to the double logarithmic data points of the normalized mean condensate size, showing an 
approximately linear behaviour. The measured slope gives the numerical estimate for the coarsening exponent Pi. To 
get sensible error estimates we slice the coarsening time window into four smaller time windows (which may overlap), 
and measure the exponent in each of the windows. The error (3i is then taken as the standard deviation of these 
measurements . 

The simulation results are shown in Figure|31for phases A and B, and in Figure0]for phase C. We plot the data on 
a double logarithmic scale for the three largest system sizes and compare to the expected scaling law given by straight 
lines. Finite size scaling of the measured exponents is given in insets. In phase C the measurements show rather large 
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FIG. 3: Verification of the predicted scaling laws, denoted by straight lines in a double logarithmic plot of the normalized 
mean condensate size for asymmetric hopping as a function of time. The predicted scaling exponents are given in 

parentheses. The insets show the finite size scaling for the numerical estimates, where filled and unfilled symbols correspond 
to the data. 

Top: Phase A with two sets of parameters. Symbols L = 512(0), 1024(A), 2048(D) for 7 = 1 and filled symbols for 7 = 0.5. 
Bottom: Phase B with two sets of parameters. Symbols L = 512(0), 1024(A), 2048(D) for 7 = 2 and L = 1024(i), 2048(a), 
4096(B) for 7 = 0.7 
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FIG. 4: Verification of the predicted scaling laws in phase C. The plots are analogous to Figure [31 except that each one 
corresponds to only one value of 7 and condensate sizes of both species are shown in filled and unfilled symbols. 
Top: Phase C with 7 = 0.5. Symbols L = 1024(0), 2048(A), 4096(D) for species 2 and filled symbols for species 1. 
Bottom: Phase C with 7 = 1.5. Symbols L = 512(0), 1024(A), 2048(D) for species 2 and filled symbols for species 1. 



errors but are in good agreement with the predictions. For 7 — 1.5 we see that, as expected, the species 1 particles 
coarsen first, but with the same exponent as species 2 particles. In phases A and B error bars are smaller, but there 
are stronger finite size effects affecting the quality of the data collapse. Nevertheless the measured scaling exponents 
are in good agreement with the predictions. 

Finite size effects strongly depend on the parameter 7, and since there are many competing mechanisms the value 
or even the sign of the resulting finite size correction is very hard to estimate. We provide a discussion of these finite 
size effects and their influence on the data collapse in the appendix. 
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Phase 


Coarsening exponents 


A 
B 

C, 7< 1 
C, 7 >1 


01 = 1/3 
Bo — 1+7 

l^ 2 2+37 
01=02=^ 
01=02 = 1/3 



TABLE II: Coarsening exponents for symmetric hopping. 




0.05 0.1 0.2 

t/T C 



FIG. 5: Verification of the predicted scaling laws for symmetric hopping in phase C with 7 = 0.5. Details of the plot are given 
in the caption to Figure [3] 

Symbols L = 512(0), 1024(A) for species 2 and filled symbols for species 1. 



C. Discussion 

One can also obtain predictions for the coarsening exponents when the hopping is symmetric. In this case, there is 
a high probability that a particle leaving a cluster site will return to the same site. The probability that it reaches 
the next cluster site is inversely proportional to the distance between cluster sites ^3 so it i s °f order L~ x . Therefore 
only every (D(L)-th excess particle will reach the next cluster site. Hence the coarsening time scales are increased by 
a factor of order L. The assumption that excess particles move independently through the bulk remains a good one 
however: the time it takes particles to enter the bulk increases by a factor O(L) (compared with the driven case) 
since most particles return to the site they have just left; this cancels the O(L) increase in the time particles spend in 
the bulk due to the diffusive rather than driven motion. Then the arguments presented for the driven case with the 
extra factor O(L) in the coarsening time scales leads to the exponents given in Table [H] We compare the prediction 
to preliminary simulation data in Figure^ where we get good data collapse on the symmetric time scale. The system 
sizes are, however, too small for a reasonable numeric estimate of the scaling exponent, but at least one can see that 
the exponent is significantly smaller than for totally asymmetric hopping (cf. Figure 0] top). 

For partially asymmetric hopping excess particles return to the cluster site they just left with non-zero probability, 
but also the probability of reaching the next cluster site in direction of the drive is of order 1, even in the limit L — > 00. 
So in contrast to symmetric hopping the coarsening time scale is only corrected by a constant factor independent of 
-L, and the coarsening exponent is the same as in the totally asymmetric case. 
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It is interesting to compare our prediction for the exponents in phase B with the results of |18| and |19( in which the 
authors study a single-species zero-range process where the hop rates, W\, . . . , wl, are site-dependent but independent 
of the particle occupation number at the departure site. They consider symmetric [l^ and asymmetric jl^ dynamics, 
where the (quenched) hop rates are drawn independently from a distribution p(w) which can be written in the form 

p(w) = [( 7 "l + 1)/(1 - a)i~ 1+1 }{w - a)^ 1 , (28) 

where w £ [a, 1] with 7, a > 0. This model undergoes a condensation transition above a critical particle density from 
a homogeneous phase to a phase with a condensate which resides at the site with the smallest hop rate. In both 
asymmetric and symmetric cases, they obtain coarsening exponents identical to those we obtain for the coarsening 
of the species 2 particles in phase B. One can think of the dynamics defined in (J3J) as a model of particles (species 2 
particles) moving on an evolving disordered background (given by the species 1 particles) . By the time the coarsening 
regime has been reached, at the cluster sites the evolving disorder is effectively quenched. Therefore it is not necessarily 
surprising that the two models exhibit similar coarsening behaviour for some distribution p(w). The reason the form 
(|28H is the relevant one for the rates JSJ is as follows. In the disordered model the coarsening is governed by the 
exchange of particles between the two slowest sites in the system. The rate at which particles are transferred between 
these two slowest sites is given by the difference between the two rates at these sites, Aw. For the distribution l|28|l. 
Aw ~ _L _ t/( 1 +t). This rate separation contributes the same factor to the coarsening time scale as that in the two- 
species model due to the dependence of the hop rate of species 2 particles on the background of species 1 particles (see 
equation 122(1 '). The remaining contributions to the coarsening time scale are then determined by the symmetry of the 
hopping dynamics, i.e. the coarsening time scale is given by a factor of order L for asymmetric dynamics, or a factor 
of order L 2 for symmetric dynamics, multiplied by the inverse rate separation. This leads to the same exponents as 
those obtained for phase B. 

V. CONCLUSION 

We have considered a two-species zero-range process which undergoes a variety of transitions to different condensate 
phases. The combination of two conservation laws and the coupling in the dynamics between the two species of 
particles leads to coarsening dynamics which are very rich compared to the single species model. In particular, we 
have considered a case in which the dynamics of one of the particle species (the species 2 particles) depends only on 
the number of particles of the other species (species 1) at the departure site, and decays to a constant value as a power 
law with exponent 7. While the stationary phase diagram discussed in Section II. B depends also on other system 
parameters, the coarsening exponents only depend (continuously) on 7 and differ from phase to phase. Further, as 
expected, the exponents depend on the symmetry of the hopping dynamics. 

Acknowledgments The authors thank the Max Planck Institute for Complex Systems, Dresden, where this work 
was initiated, for hospitality. TH was supported by EPSRC programme grant GR/S10377/01. 

APPENDIX 
DISCUSSION OF FINITE SIZE EFFECTS 

In the following we discuss qualitatively the most basic mechanisms leading to finite size effects and illustrate how 
they depend on the system parameters. We start with two competing effects which influence the bulk density of 
the condensing species, connected with the definition of cluster sites by a threshold, (i) Excess particles exchanged 
between cluster sites increase the bulk density in finite systems. In phase B, for example, this effect is of order 
£1-7/(1+7) = by reasoning given in Section III.B and the expansion l|22l) . It decreases with increasing 7 

and leads to a decrease of the mean condensate size. This effect is shown for 7 < 1 in Figure HO (bottom left), where 
P2,buik is plotted, decreasing towards pi, c with increasing system size L. (ii) On the other hand, bulk fluctuations 
increase with 7 and in finite systems they can exceed the threshold for cluster sites. This leads to an increase of 
the number of condensed particles, or a decrease of /02,buik, which dominates over effect (i) for 7 > 1. This can be 
seen in Figure (bottom right) where p2,buik now increases towards p2, c - Since these fluctuations are typically small 
compared to other cluster sites, this leads also to a decrease of the mean condensate size, (iii) Further, in finite 
systems the nucleation and the coarsening regime are not clearly separated but overlap to a large extent, as can also 
be seen in Figure (bottom). Ongoing nucleation effectively slows down the coarsening, leading to a decrease of the 
coarsening exponents for finite system sizes, as can be seen in the insets of Figures 3 (bottom) and 4 (top). 
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FIG. 6: Finite size effects in phase B. /9i,buik shown on the top is smaller than its limiting value pi = 0.5 for L — » oo since on 

species 2 cluster sites kx ~ L 1/(1+7) . With increasing 7 this effect becomes weaker. pa.bulk shown on the bottom converges to 

its limiting value p2, c from above (for 7 < 1) and from below (for 7 > 1) as explained in the text. 

Left: 7 = 0.7, 6 = 1, c = 3, pi = 0.5, p2 = 5, critical density p2, c = 1-28 

Symbols: L = 256(i), 512(*), 1024(B), 2048(a), 4096(«) 

Right: 7 = 2, b — 1, c = 4, pi = 0.5, P2 = 10, critical density p2, c = 3.40 

Symbols: L = 256(*), 512(*), 1024(B), 2048(a) 



On top of these corrections, which are also present in single species models, there are mechanisms specific to 
two-component systems, (iv) The coarsening dynamics does not only depend on the occupation number of the 
condensing species, but on the relation between ki and &2- As discussed above, condensates with different relations 
have different lifetime. So in the limit L — > 00 ratios differing by some factor L a are dynamically separated and 
only one relation dominates the coarsening. But for finite systems, ratios with shorter lifetimes also contribute to 
the observed behaviour. Thus data for single species systems, where this effect does not occur, are generically better 
than our data (cf. a 2j). (v) Finally, we consider a phenomenon specific for phase B. As discussed above, on species 2 
cluster sites there are of order species 1 particles due to compatibility with the bulk. For 7 < 1 this is larger 

than a typical bulk fluctuation of order L 1 / 2 , and reduces the bulk density of species 1 particles, whereas for 7 > 1 
the effect is much weaker. This can be seen in Figure [S] (top) , where pi^uik is plotted. 

The variety of effects leads to a diverse behaviour and a quantitative prediction of the finite size corrections seems 
to be not feasible. However, it is straightforward to numerically fit the leading order corrections for the prefactor and 
the exponent of the scaling law @ , using the ansatz 

( m *(t)) L /~r /1 1 /-» ,rS^(t\^ +C3/L 



{pi - Pi,c)L 



C 1 {l + C 2 /L^)(-) . (29) 



Consider for instance the data in phase C with 7 = 0.5 given in Figure 0] (top). The finite size corrections in the inset 
suggest 5^ = 1 and the best fit values for the other parameters are 8\ = 0.69, C\ — 0.51, C2 = 2.3, C3 = —21 for 
species 1 and 81 = 0.60, C\ = 0.59, C2 = 2.1, C3 = —13 for species 2. In FigureQwe plot the corrected data 



( mi (t)) f L s (m(t)) L ,t\-c 3 /L s * 



{pi - Pi,c)L {pi - Pi, c )L 



(d + Ca/L 51 ) 1 , (30) 




and the data collapse improves drastically. 
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